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Hard X-ray and extremely low energy bulk-sensitive photoelectron spectroscopy has been 
performed in the temperature range of 100-330 K for FeaCU. In the high temperature phase 
just above the Verwey transition, the intensity at the Fermi level (Ef) is still negligible, but 
it increases gradually with further increasing the temperature (250 K, 330 K) in consistence 
with the temperature dependence of the conductivity. The spectral behaviors near Ef with 
temperature are well explained by the model, which takes the polaron effect into account. 

KEYWORDS: strongly correlated electron system, Verwey transition, Magnetite, 
photoelectron spectroscopy, polaron 



Fe3C>4 (magnetite) is known as one of the oldest mag- 
netic oxides with high potential for applications to spin- 
electronics (spintronics). Many theoretical results have 
predicted the half metallic ferromagnetism (HMF) in 
Fea04 with a conductivity by minority spin electrons and 
the semiconducting behavior by the majority spin chan- 
nel. 1-4 ) Recently, the experimental evidence for the HMF 
is given for the epitaxial Fe304 (111) thin films at room 
temperature (RT) by means of the spin-resolved photo- 
electron spectroscopy in the photon energy (his) range of 
20-60 eV. 5 ' 6 ) 

Fe304 is well known to show the Verwey transition 
across Ty ~ 123 K, 7 ) where the first order transition of 
conductivity takes place with its decrease by about two 
orders of magnitude accompanied with the lattice distor- 
tion. 8 ) There are controversial discussions on the origin 
of this transition. 9 ~ 14 ) For example, Fe3C>4 is thought to 
be metallic above Ty, although the temperature depen- 
dence of the conductivity in Fe3C>4 suggests the insulat- 
ing character below ^300 K. The photoelectron spectral 
weight near the Fermi level (Ep) has been found to be 
very weak even above Ty 15 ~ 17 ' in contrast to the results 
of the band calculation. The insulator to insulator tran- 
sition across Ty with increasing the temperature is pro- 
posed by conventional photoelectron spectroscopy (PES) 
performed in the hv range of 20-100 eV. 15 > On the other 
hand, the insulator to metal transition (IMT) is proposed 
by a high resolution photoemission study by Chainani et 
al. 16 ' Very recently more bulk sensitive soft X-ray PES 
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(SXPES) has been performed at hv = 707.6 eV under 
the Fe 2p-3d resonance excitation, and revealed the in- 
sulator to insulator transition across Ty. 17 ) 

As well known nowadays, the surface electronic struc- 
tures are much different from the bulk electronic struc- 
tures 18 ~ 20 ) in strongly correlated electron systems be- 
cause U/t (U : electron correlation energy, t : intersite 
hopping energy) is much larger in the surface region com- 
pared with that in the bulk region. Even in the SX- 
PES, however, the surface spectral weight near Ep can 
be as high as 20-40 %. Therefore more bulk-sensitive 
PES studies are desired for a long time. In this decade, 
hard X-ray PES (HAXPES) has widely been recognized 
as a highly bulk-sensitive technique because of the long 
inelastic mean-free path (IMFP) and intensively applied 
to various strongly correlated electron systems. 21 ~ 23 ) The 
IMFP reaches up to ~100 A at hv ~ 8 keV. 24) Meanwhile 
extremely low energy PES (ELEPES) with a laser excita- 
tion near 7 eV turned out to provide rather bulk sensitive 
results. 25 ' 26 ) The IMFP for ELEPES is expected to reach 
values comparable to that for HAXPES under certain 
conditions. We have here employed both HAXPES and 
ELEPES for studying the electronic structures of Fe3C>4 
across the Verwey transition and in high temperatures 
up to 330K, where we confirmed that the behavior of 
bulk PES can be consistently understood by considering 
the polaron effects. 

The experiments were performed on single crystalline 
Fe3 04 grown by the floating-zone method. HAXPES 
measurements were performed at BL19LXU in SPring- 
8 27) with an MB Scientific (MBS) Al-HE hemispherical 
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analyzer. The ELEPES measurements were performed 
by use of synchrotron radiation at BL7U in UVSOR 28 ) 
with an MBS-A1 analyzer. The photon energy employed 
in HAXPES (ELEPES) was 7941 eV (7.5 eV) and the en- 
ergy resolution was set to ^120 meV (<10 meV). Clean 
surfaces were obtained by fracturing in situ at 140 K 
(room temperature) in HAXPES (ELEPES). The Fermi 
level was calibrated by the Fermi edge of Au electrically 
connected to the sample. The HAXPES (ELEPES) mea- 
surements were performed at 100, 140, 250, and 330 K 
(140, 250, and 330 K). 

Figure 1(a) displays the temperature dependence of 
the HAXPES valence-band spectra of Fe304 in both low 
temperature (LT) phase (100K) and high temperature 
(HT) phase (140, 250, and 330 K) with respect to T v 
measured with the resolution of ^200 meV. The broad 
band ranging from 5 to 12 eV is ascribed to the so-called 
O 2p valence band. One notices that the intensity near 
8 eV is much stronger than the intensity near 4 eV in 
strong contrast to the results of low energy and soft X- 
ray PES. 15-17 ) According to the photoionization cross 
section 29 ^ the strong intensity near 8 eV is ascribed to the 
Fe 4s component hybridized with the O 2p component 
similar to vanadium oxides, 30, 31 ) whereas the structures 
in the binding energy (Eb) below 5 eV are either Fe 3d- 
O 2p anti-bonding state or the O 2p non-bonding state. 
It should be noticed that the transition metal 4s + O 
2p components in 5-12 eV in the HAXPES spectra have 
less temperature dependence than those for most vana- 
dium oxides with metal-to-insulator transitions. 30, 31 ) As 
shown in Fig. 1(b), however, a clear shift of the spectral 
weight is observed near the Fermi level (Ep) across the 
Verwey transition in the spectra measured with higher 
resolution of 120 meV. In the LT phase (100 K), the 
spectrum has the small energy gap as judged from the 
position of the \ow-Eb threshold and negligible intensity 
at Ep, as expected for an insulator. In the HT phase just 
above Ty (at 140 K), the spectral onset becomes shifted 
towards Ep, but the intensity just at Ep is still very 
weak as already seen in SXPES. 17 - 1 Furthermore, spectral 
change was observed with further increasing temperature 
up to 330 K in HT phase. The integrated spectral inten- 
sity between -0.1 and 0.1 eV is plotted as a function of 
temperature in Fig. 1(c) by the blue triangle marks. It 
increases gradually with temperature in HT phase and 
the enhancement of the spectral weight between 250 and 
330 K is smaller than that between 140 K and 250 K. 
For a more quantitative discussion, we phenomenolog- 
ically define the spectral onset E on as the intersection 
of the leading edge with the zero intensity base line as 
shown in the inset of Fig. 1(c). E on is located below Ep 
at 100 K in the LT phase and just above Ty (at 140K) 
as expected for an insulator, whereas it is above Ep at 
250 K and 330 K. 

In order to see the detailed temperature dependence of 
the intensity at Ep in the HT phase, ELEPES has been 
performed in the temperature range of 140-330 K for 
Fea04 with much better energy resolution (<10 meV). 
Figure 2 displays the temperature dependence of the 
ELEPES spectra near E F of Fe 3 4 in the HT phase. The 
spectra are normalized by the area with Eb larger than 



Fig. 1. (color online) (a) Temperature dependence of the hard 
X-ray photoemission (HAXPES) spectra of the valence band of 
Fe3C>4 in the wide Eb region normalized by the area above Eb = 
3 eV. (b) Temperature dependence of the HAXPES valence band 
spectra of Fe304 near Ep measured with the total resolution of 
120 meV. The dots and solid lines show the experimental and 
the smoothed data, respectively, (c) Temperature dependence of 
the energy (E on ) of the spectral onset (red diamonds) (see inset 
for definition) and the spectral intensity at Ep (blue triangles) 
above Ty tentatively obtained by integrating the photoelectron 
intensity from -0.1 to 0.1 eV. 

1.0 eV since this region dominated by the O 2p states 
had no temperature dependence in HAXPES spectra. 
The shapes of ELEPES spectra are much different from 
those of HAXPES spectra. For example, the intensity 
increases still above 0.8 eV in contrast to the HAXPES. 
According to the photoionization cross section, the ratio 
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Fig. 2. (color online) Temperature dependence of the ELEPES 
spectra for Fe304 in the high-temperature phase. The inset 
shows the expanded ELEPES spectra near Ep. 

between the O 2p and Fe 3d components is much larger in 
ELEPES spectra than in the HAXPES. 29 ) Therefore, the 
spectral shapes of ELEPES mainly reflect the O 2p com- 
ponent hybridized with the Fe 3d states. At 140 K, the 
intensity is very weak at Ep, but it becomes finite and in- 
creases with temperature (inset of Fig. 2) in consistence 
with the HAXPES result. The temperature dependence 
of the intensity at Ef can not be simply explained by the 
Femi-Dirac distribution function. Moreover, the intensity 
at around 0.6 eV is gradually suppressed with increasing 
temperature. These temperature dependence clearly in- 
dicates that the electronic states of FC3O4 change even 
within the HT phase. These spectral changes can well 
explain the unusual thermal behavior of the conductiv- 
ity above Ty, which increases toward 330 K. One notices 
here that the sharp Fermi edge is never observed in FC3O4 
at 330K. In the LDA calculation, 3 ) however, large den- 
sity of states (DOS) is predicted near Ef due to the Fe 
3d (hybridized with O 2p) down-spin states in contrast 
to the experimental results. 

This discrepancy between the PES spectra and LDA 
calculation results can be explained by the formation of 
the polaronic quasiparticle due to the electron-phonon 
coupling, which is often discussed in HT phase of Fe 3 4 . 
Degiorgi et al. explained the unusual temperature be- 
havior of the conductivity above Ty by using the small 
polaron model which also takes the polaronic short range 
order into account. 32 ) The mid-infrared polaron peak 
predicted from a polaronic binding energy of e p ~ 300 
meV has really been observed in the optical conductivity 
for Fe3 04. 33 ' 34 ) Strong coupling to the lattice vibration 
and the formation of small polarons have also been in- 
ferred from the small magnitude of carrier mobility. 35 ^ 
Besides, the temperature dependence of the HAXPES 
and ELEPES spectra near Ep of Fe3 04 is very similar to 
the spectral changes of (TaSe4)2l and K0.3M0O3, which 
are known as the quasi-one-dimensional Peierls systems 
with polaronic effects. 36 ' 37 ^ The previous SXPES mea- 
surements have also suggested the formation of the small 



Fig. 3. (color online) (a) Raw HAXPES spectrum near Ep (blue 
curve) and the background-subtracted spectrum (red curve) of 
Fe 3 4 at 140 K (above Ty). (b) Model fit of the background- 
subtracted HAXPES spectra (see text for the label A and B). 
(c) Model fit of the ELEPES spectra with the same parameters 
as HAXPES spectra. Inset shows the temperature dependence 
of the dimensionless electron-phonon coupling constant g 2 . 

polaron in HT phase of Fes04 by using of the theoreti- 
cal polaronic model of Alexandrov and Ranninger, 38 ) in 
which the electron removal spectrum of a system with 
strong electron phonon coupling can be written as 

00 2n 

I(E B ) oze-3 2 N p (E B ) + Y / e~ 92 ^N p (E B -nLU ). (1) 

The first term describes the polaronic quasiparticle 
band (N p (Eb)) renormalized by a factor e~ 9 , where 
g 2 = s p /ujo is a dimensionless electron-phonon coupling 
strength given by e p (polaronic binding energy) divided 
by ujq (characteristic phonon energy). The spectral in- 
tensity is partially transferred to multiple phonon side 
bands shifted by nuj toward higher binding energies. If 
the coupling constant g 2 becomes stronger, the polaron 
quasiparticle weight is exponentially suppressed in the 
photoemission spectrum. This polaronic model could fit 
the experimental spectra of SXPES at 180 K with pa- 
rameters of g 2 ~ 5 and uj ~ 70 meV in consistence with 
both the polaronic binding energy (^300 meV) and the 
effective mass (~100~200mo) derived from the infrared 
spectroscopy. 17 ) 

We have tried to analyze the HAXPES and ELEPES 
results by modeling the above polaronic spectral function 
based on eq. (1) and the LDA calculation. 3 ) The LDA 
DOS near Ep consists of two components, one cross- 
ing Ep composed of minority spin states and the other 
one centered at Eb — 0.6 cV composed of majority spin 
states. The polaronic model is applied only to the for- 
mer component, because the spectral change around Eb 
= 0.6 eV is much smaller than that near Ep. In order 
to reproduce the HAXPES (ELEPES) spectra in Fe 3 4 , 
we have taken the photoionization cross section ratio be- 
tween Fe 3d and O 2p states into account 39 ) and broad- 
ened the theoretical spectrum by 120 meV (10 meV) to 
consider the instrumental resolution. We have first sub- 
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tracted the Fe 4s + 2p contribution above Eb = 1 cV 
in HAXPES as shown in Fig. 3(a) by the dashed curve 
and then approximated the component centered at Eb 
= 0.6 eV by a proper Gaussian peak. 40 ) According to 
the results of the infrared study for Fe3C>4 two phonon 
modes were reported, 33, 34 ' one at ^40 meV and the other 
at ^70 meV. We have tried to fit the spectra with either 
phonon energy because it is not known which phonon 
mode couples to the electron. A good agreement between 
the calculation and the experiments is obtained by using 
the higher phonon energy of 70 meV as in the case of SX- 
PES results. Therefore, the characteristic phonon energy 
ujq is fixed to 70 meV. The coupling constant g 2 is treated 
as the fitting parameter. We note that the temperature 
dependence of the HAXPES and ELEPES spectra in HT 
phase could not be satisfactory explained without chang- 
ing g 2 with temperature as discussed below. 

Figure 3 (b) and (c) show the fitting results of 
the background-subtracted HAXPES spectra and the 
ELEPES spectra. In Fig. 3 (b), the dashed curves (Peak 
A) show the spectra derived through eq. (1) from the 
DOS component crossing Ep and the dot-dashed curves 
(Peak B) show the spectra from the DOS centered at Eb 
= 0.6 eV. We can well reproduce the spectra near Ef at 
each temperature with the common coupling strength 
g 2 of 3.4, 2.2 and 1.8 at 140, 250 and 330 K through 
HAXPES and ELEPES. Such a temperature-dependent 
coupling strength has been implied by the optical con- 
ductivity. 33, 34 ' The coupling strength g 2 is plotted as a 
function of temperature in the inset of Fig. 3 (c). At 
140 K, e~ 3 - 4 ~ 0.033 means the suppression of the spec- 
tral weight at Ep of more than 96 %, where the value 
of the g 2 — 3.4 [e p — 238 meV being consistent with 
the polaronic binding energy (^300 meV) derived from 
the infrared spectroscopy 33,34 '] reflects the almost insu- 
lating feature and the formation of the small polaron. 
Then g 2 (e p ) gradually decreases with increasing tem- 
perature. The suppression of the spectral weight due to 
the formation of the polaron becomes reduced as e~ 2,2 
- 0.11 and e- 1 - 8 ~ 0.17 around 250 -330 K, indicating 
the variation from the small polaron to the large polaron 
with increasing temperature. It has been suggested from 
the results of neutron scattering, 41 ' /zSR, 42 ' transport 
phenomena, 43 ' and Mossbauer spectroscopy 44 ' that the 
change of the electronic states takes place below around 
250-300 K corresponding to the formation of the large 
polaron. Our results strongly support the scenario as the 
variation takes place from the small polaron to the large 
polaron with increasing temperature in the HT phase of 
Fe 3 4 . 

In conclusion, we have performed HAXPES and 
ELEPES for Fea04 and observed the detailed spectral 
changes with increasing the temperature. In the HT 
phase just above the Verwey transition, the intensity at 
Ep is still negligible. However, it increases gradually with 
increasing temperature. This temperature dependence 
can be well explained by considering the polaron effect. 
Since the polaronic coupling strength becomes weaker at 
high temperatures far from Ty, the intensity near Ep in- 
creases inducing the increased conductivity of Fea04 at 
higher temperatures. 
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